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SUBSTITUTE SPECIFICATION 



Attorney Docket No. 36856.1330 

ELECTROMAGNETIC FIELD ANALYZER, ELECTROMAGNETIC FIELD 
ANALYZING PROGRAM, AND STORAGE MEDIUM FOR RECORDING THE 

PROGRAM 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to electromagnetic field 
simulation technology used in developing and designing high- 
frequency products such as multilayer products, giga filters, 
and EMI (Electro-Magnetic Interference) elimination filters, 
and more particularly, the present invention relates to an 
electromagnetic field analyzer, an electromagnetic field 
analyzing program, and a storage medium for recording the 
program in which the simulation time is reduced. 

2. Description of the Related Art 

In recent years, the development and design of high- 
frequency products have been aggressively pursued, and, in 
order to perform an electromagnetic field analysis of high- 
frequency products, electromagnetic field simulation software 
has been commonly used. Generally, there are many instances 
where a finite element method is used for electromagnetic 
field analysis. In the electromagnetic field simulation using 
the finite element method, since it takes a long time to solve 
simultaneous linear equations, high-speed solving methods are 
desired. In particular, in the field of electromagnetic field 
analysis, since only a direct method called a Gauss' 

elimination method can be used, the problem of long 
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calculation times is more serious. 

However, a paper (R. Hiptmair, "Multigrid Method for 
Maxwell's Equations, " SIAM Journal of Numerical Analysis, vol. 
36, no. 1, pp. 204-225, 1999), which is referred to herein as 
Non-patent document 1 , and which was made public by Hiptmair 
in 1999, proved that a solution method called a multi-grid 
method can be also used in an electromagnetic field analysis 
and the method greatly reduces the calculation time as 
compared to a direct method. 

In the multi-grid method, an analysis object is divided 
into two classes of fine elements and coarse elements, and the 
solution of fine elements is determined using the solution of 
coarse elements. 

When simultaneous linear equations are calculated by the 
Gauss' elimination method, generally the calculation time is 
proportional to the third power of the dimensions of a matrix. 
Since the size of the dimensions corresponds to the number of 
elements, when the number of elements is doubled, the 
calculation time is 8 times as large, and, when the number of 
elements is 10 times as large, the calculation time is 1,000 
times as large. When the Gauss' direct method is applied to 
only coarse elements, an approximate solution to fine elements 
is obtained by using the solution to coarse elements, and the 
Gauss' elimination method is not used, the calculation time 
for fine elements can be reduced to 1/8 and 1/1000, 
respectively. Since additional calculation time is added, 
although the practical calculation time is not reduced to such 
a great extent, the calculation time is always greatly reduced. 

The multi-grid method can be broken down into two methods: 
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a method using a nested mesh and a method using a non-nested 
mesh. However, since a prolongation matrix is incomplete as 
described later, the method using a nested mesh must be used. 

Figs. 1 2A and 1 2B show the division into elements using a 
nested mesh. Fig. 1 2A shows the division into elements using 
a nested mesh in the case of two-dimensional analysis where a 
coarse triangular element is uniformly divided to form fine 
triangular elements. Furthermore, Fig. 1 2B shows the division 
into elements using a nested mesh in the case of three- 
dimensional analysis where a coarse tetrahedral element is 
uniformly divided to form eight fine tetrahedral elements. As 
understood from Figs. 1 2A and 12B, there are geometrical 
restrictions between coarse elements and fine elements when 
using a nested mesh. 

Another prior art method is described in non-patent 
document 2: D. Dibben and T. Yamada, "Non-Nested Multigrid and 
Automatic Mesh Coarsening for High Frequency Electromagnetic 
Problems" IEEJ Investigating Research Committee Material, SA- 
02-34, pp. 71-75, 2002. 

Figs. 1 3A and 1 3B describe a first problem that occurs 
when using a nested mesh. Fig. 1 3A shows how a circle is 
divided into elements using a nested mesh to form coarse 
elements. Furthermore, Fig. 1 3B shows how coarse elements are 
divided into fine elements using a nested mesh. As understood 
from Figs. 1 3A and 13B, a curved surface cannot be correctly 
expressed because of the geometrical restrictions of a nested 
mesh . 

Furthermore, Fig. 1 3C shows how a circle is divided into 
elements using a non-nested mesh to form fine elements. As 
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understood from Fig. 13C, since there are no geometrical 
restrictions when using a non-nested mesh, a curved surface 
can be correctly expressed. 

Figs. 1 4A to 14C describe a second problem that occurs 
when using a nested mesh. Fig. 1 4A shows how a square is 
divided into elements using a nested mesh to form coarse 
elements. Furthermore, Fig. 1 4B shows how coarse elements are 
divided into fine elements. As understood from Figs. 1 4A and 
14B, coarse elements cannot be partially divided into fine 
elements because of the geometrical restrictions of a nested 
mesh . 

Furthermore, Fig. 14C shows how a square is divided into 
elements using a non-nested mesh to form fine elements. As 
understood from Fig. 14C, coarse elements can be partially 
divided into fine elements because there are no geometrical 
restrictions when using a non-nested mesh. 

Most industrial products have curved surfaces, such as a 
circular cylinder and square surfaces, and they inevitably 
contain partially fine portions. Accordingly, it is difficult 
to apply a nested mesh to industrial products. On the other 
hand, since there are no geometrical restrictions in a non- 
nested mesh, fine elements can be freely formed. However, 
since the positional relationship is not systematic, it is 
difficult to make the electromagnetic field of coarse elements 
related to the electromagnetic field of fine elements. It is 
a prolongation matrix to make coarse elements related to fine 
elements, but no precise prolongation matrix has been found. 
This means that a multi-grid method cannot be practically used, 
and accordingly, nothing can be used except for a Gauss' 
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direct method, which requires a lot of time, in the analysis 
of magnetic field. 

SUMMARY OF THE INVENTION 
To overcome the problems described above, preferred 
embodiments of the present invention provide an 
electromagnetic field analyzer, an electromagnetic field 
analyzing program, and a recording medium having the program 
recorded thereon, in which the analysis of an electromagnetic 
field can be performed via a multi-grid method using a non- 
nested mesh, and an electromagnetic field analyzer, an 
electromagnetic field analyzing program, and a recording 
medium in which the program has been recorded, in which the 
time required for an electromagnetic field analysis is greatly 
reduced. 

According to a preferred embodiment of the present 
invention, an electromagnetic field analyzer includes dividing 
means for dividing form data as an analysis object into coarse 
elements and fine elements, forming means for forming a matrix 
defining an electromagnetic field vector of the coarse 
elements divided by the dividing means related to an 
electromagnetic field vector of the fine elements, and 
calculating means for calculating an approximate solution of 
the electromagnetic field vector of the fine elements by 
applying an iteration method of simultaneous linear equations 
while referring to the matrix formed by the forming means. 

Preferably, the forming means forms the matrix by 
expressing the elements of electromagnetic field vectors at 
the sides at a fine element using an interpolation function in 
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the coarse elements. 

More preferably, when the length of a side li of the fine 
element is |ljj, an interpolation function showing the 
relationship between the electromagnetic field at position x 
of the side li of the fine element and the electromagnetic 
field at a side j of the coarse element is Nj C (x) , and the 
tangential vector of the side li of the fine element is ti, the 
forming means forms a matrix Pij preferably by Expression (14) 
to be described later. 

More preferably, in the calculating means, the accuracy of 
an approximate solution of the electromagnetic field vector of 
the fine elements is improved such that high-frequency 
components included in the approximate solution of the 
electromagnetic field vector of the fine element are removed 
by applying a stationary iteration method of simultaneous 
linear equations, such that a residual in the fine elements is 
mapped to a residual in the coarse elements by using the 
matrix formed by the forming means, such that a correction 
vector to the coarse elements is formed by applying a direct 
method or a non-stationary iteration method of simultaneous 
linear equations, and such that a correction vector to the 
fine elements is obtained by using the matrix formed by the 
forming means. 

According to another preferred embodiment of the present 
invention, in a computer program to control a computer to 
execute an electromagnetic field analyzing method for 
analyzing an electromagnetic field as an analysis object, the 
computer includes a first storage device arranged to store 
divided elements and a second storage device arranged to store 
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a matrix, and the electromagnetic field analyzing method 
includes the steps of dividing form data as an analysis object 
into coarse elements and fine elements and storing the 
elements in the first storage device, forming a matrix 
defining an electromagnetic field vector of the coarse 
elements stored in the first storage device related to an 
electromagnetic field vector of the fine elements and storing 
the matrix in the second storage device, and calculating an 
approximate solution of the electromagnetic field vector of 
the fine elements by applying an iteration method of 
simultaneous linear equations while referring to the matrix 
stored in the second storage device. 

According to another preferred embodiment of the present 
invention, in a computer readable storage medium having stored 
thereon a program to control a computer to execute an 
electromagnetic field analyzing method for analyzing an 
electromagnetic field as an analyzing object, the computer 
includes a first storage device arranged to store divided 
elements and a second storage device arranged to store a 
matrix, and the electromagnetic field analyzing method 
includes the steps of dividing form data as an analysis object 
into coarse elements and fine elements and storing the coarse 
and fine elements in the first storage device, forming a 
matrix defining an electromagnetic field vector of the coarse 
elements stored in the first storage device related to an 
electromagnetic field vector of the fine elements and storing 
the matrix in the second storage device, and calculating an 
approximate solution of the electromagnetic field vector of 
the fine elements by applying an iteration method of 
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simultaneous linear equations while referring to the matrix 
stored in the second storage device. 

Since the forming means forms a matrix defining an 
electromagnetic field vector of coarse elements divided by the 
dividing means related to an electromagnetic field vector of 
fine elements, it is possible to perform an electromagnetic 
field analysis via a multi-grid method using a non-nested mesh. 
Furthermore, since the calculating means calculates an 
approximate solution of the electromagnetic field vector of 
fine elements while referring to the matrix formed by the 
forming means, the calculating means is not required to 
calculate the electromagnetic field of fine elements by using 
a direct method. Thus, the time required for an 
electromagnetic field analysis is greatly reduced. 

Furthermore, since the forming means forms a matrix such 
that an electromagnetic field vector element at a side of fine 
elements is expressed by using an interpolation function in 
coarse elements, a matrix is easily formed. 

Furthermore, in the forming means, since a matrix Pij is 
formed preferably by Expression (14) to be described later, a 
matrix is even more easily formed. 

Furthermore, in the calculating means, an approximate 
solution of an electromagnetic field vector of fine elements 
is corrected by using a matrix formed by the forming means. 
Thus, the accuracy of the approximate solution of an 
electromagnetic field vector of fine elements is greatly 
improved at high speed. 

Other features, elements, steps, characteristics and 
advantages of the present invention will become more apparent 
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from the following detailed description of preferred 
embodiments thereof with reference to the attached drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows an electromagnetic field analyzer according 
to a preferred embodiment of the present invention. 

Fig. 2 is a block diagram showing the construction of the 
electromagnetic field analyzer according to a preferred 
embodiment of the present invention. 

Fig. 3 is a block diagram showing the functional structure 
of the electromagnetic field analyzer. 

Fig. 4 is a flow chart for describing the procedure of the 
electromagnetic field analyzer. 

Fig. 5A shows a microstrip line divided into large 
elements . 

Fig. 5B shows a microstrip line which is analyzed using 
the large elements shown in Fig. 5A, and in which the elements 
are subdivided based on their error estimation. 

Fig. 6 shows components of an electromagnetic field that 
are in contact with the sides of a tetrahedron. 

Fig. 7 shows the relationship between an electromagnetic 
field vector u F in Q F and an electromagnetic field vector u c in 

a c . 

Fig. 8 shows one example of a rectangular waveguide as an 
analysis object of an electromagnetic field analyzer according 
to a preferred embodiment of the present invention. 

Fig. 9 shows the calculation time to the number of 
elements when the rectangular waveguide shown in Fig. 8 is 
driven in TE10 mode of 2.45 GHz, which is compared between the 
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case where a related direct method is used and the case where 
a method of a preferred embodiment of the present invention is 
used. 

Fig. 10 shows one example of a patch antenna as an 
analysis object of the electromagnetic field analyzer 
according to a preferred embodiment of the present invention. 

Fig. 11 shows the calculation time to the number of 
elements when the patch antenna shown in Fig. 10 is driven at 
7.0 GHz, which is compared between the case where a related 
direct method is used and the case where a method of the 
present invention is used. 

Fig. 1 2A shows elements divided by using a nested mesh in 
the case of two-dimensional analysis. 

Fig. 1 2B shows elements divided by using a nested mesh in 
the case of three-dimensional analysis. 

Fig. 1 3A shows coarse elements formed such that a circle 
is divided into elements by using a nested mesh. 

Fig. 1 3B shows fine elements into which coarse elements 
are divided by using a nested mesh. 

Fig. 1 3C shows fine elements formed such that a circle is 
divided into elements by using a non-nested mesh. 

Fig. 1 4A shows coarse elements formed such that a square 
is divided into elements using a nested mesh. 

Fig. 1 4B shows fine elements into which coarse elements 
are divided by using a nested mesh. 

Fig. 14C shows fine elements formed such that a square is 
divided into elements by using a non-nested mesh. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Fig. 1 shows an electromagnetic field analyzer according 
to a preferred embodiment of the present invention. The 
electromagnetic field analyzer preferably includes a computer 
main body 1, a display device 2, an FD (Flexible Disk) drive 3 
to which an FD 4 is attached, a keyboard 5, a mouse 6, a CD- 
ROM (Compact Disk- Read Only Memory) drive 7 to which a CD-ROM 
8 is attached, and a network communications device 9. An 
electromagnetic field analyzing program is supplied through a 
recording medium such as an FD 4, or a CD-ROM 8, or other 
suitable media. An electromagnetic field analysis is 
performed such that the electromagnetic field analyzing 
program is executed by the computer main body 1. Furthermore, 
the electromagnetic field analyzing program may be supplied to 
the computer main body 1 by another computer through the 
network communications device 9. 

Fig. 2 is a bock diagram showing the construction of the 
electromagnetic field analyzer according to a preferred 
embodiment of the present invention. The computer main body 1 
shown in Fig. 1 includes a CPU (Central Processing Unit) 10, a 
ROM (Read Only Memory) 11, a RAM (Random Access Memory) 12, 
and a hard disk 13. The CPU 10 performs processing while the 
CPU 10 inputs and outputs data to the display device 2, the FD 
drive 3, the keyboard 5, the mouse 6, the CD-ROM device 7, the 
network communications device 9, the ROM 11, the RAM 12, or 
the hard disk 13. The electromagnetic field analyzing program 
in the FD 4 or the CD-ROM 8 is preferably stored on the hard 
disk 13 through the FD drive 3 or the CD-ROM drive 7. The 
electromagnetic field analysis is performed such that the CPU 
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10 properly loads the electromagnetic field analyzing program 
from the hard disk 13 in the RAM 12 and executes the program. 

Fig. 3 is a block diagram showing the functional structure 
of the electromagnetic field analyzer according to a preferred 
embodiment of the present invention. The electromagnetic 
field analyzer includes a form data storage portion 21 for 
storing form data as an analysis object, a division-into- 
element portion 22 for dividing the form data, such that the 
size of elements is changed, an element storage portion 23 for 
storing the elements divided by the division-into-element 
portion 22, a vector /matrix forming portion 24 for forming a 
vector and matrix in simultaneous linear equations, a 
prolongation matrix forming portion 25 for forming a 
prolongation matrix, a vector/matrix storage portion 26 for 
storing the vector and matrix formed by the vector /matrix 
forming portion 24, a prolongation matrix storage portion 27 
for storing the prolongation matrix formed by the prolongation 
matrix forming portion 25, an approximate value calculation 
portion 28 for calculating an approximate value of an 
electromagnetic field vector, an approximate solution 
correction portion 29 for correcting the approximate solution 
calculated by the approximate value calculation portion 28, 
and an electromagnet ic- field vector storage portion 30 for 
storing, as an electromagnetic field vector, the approximate 
solution after having been corrected by the approximate 
solution correction portion 29. 

Moreover, the form data storage portion 21, the element 
storage portion 23, the vector /matrix storage portion 26, the 
prolongation matrix storage portion 27, and the 
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electromagnetic field vector storage portion 30 are provided 
in a fixed area in the RAM 1 2 or on the hard disk 1 3 shown in 
Fig. 2. Furthermore, the functions of the division- into- 
element portion 22, the vector /mat rix forming portion 24, the 
prolongation matrix forming portion 25, the approximate value 
calculation portion 28, and the approximate solution 
correction portion 29 are arranged such that the CPU 10 shown 
in Fig. 2 executes the electromagnetic field analyzing program 
loaded into the RAM 12. 

Fig. 4 is a flow chart for describing the procedure of the 
electromagnetic field analyzer according to a preferred 
embodiment of the present invention. First, the division- 
into-element portion 22 reads form data as an analysis object 
stored in the form data storage portion 21 , divides the form 
data into elements, and stores the divided elements in the 
element storage portion 23 (S11). The division-into-element 
portion 22 divides the form data into triangular or 
quadrangular elements in the case of two-dimensional analysis, 
and divides the form data into tetrahedral, triangle-pole- 
shaped, or square-pole-shaped elements in the case of three- 
dimensional analysis . 

The division-into-element portion 22 divides the form data 
while changing the size of the elements. In order to make the 
description simple, elements of two different sizes are formed 
and it is assumed that a set of larger (coarse) elements is 
made Q c and a set of smaller (fine) elements is made Q F . 
Moreover, a superscript of c represents coarse, and a 
superscription of F represents fine. 

The sets of Q. c and Q F may be independently formed, however, 
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when the set of Q. c is first formed and then the set of Q F is 
formed by subdividing the elements, the processing time is 
reduced. Regarding the method of subdivision, there is a 
method for subdividing elements which are larger in volume and 
a method for analyzing elements by using £l c and subdividing 
elements having a large error. Either of them may be used, 
however, the method in which error is evaluated has better 
calculation accuracy. Moreover, regarding the method for 
evaluating error, Reference Document No. 1 (Norio Kamiya et 
al., "Computer and Application's Monk, Special Edition/Error 
Estimation in Software and Adaptive Element," No. 42, CORONA 
PUBLISHING CO., 1993) may be referred to. 

Figs. 5A and 5B show examples of division into elements in 
a microstrip line. In Fig. 5A, the microstrip line is divided 
into larger elements. Furthermore, in Fig. 5B, an analysis is 
performed by using the larger elements in Fig. 5A and the 
elements are subdivided based on their error estimation. 
Since the electromagnetic field is concentrated at the edge of 
an electrode, the portion is subdivided into smaller elements. 

Next, the vector/matrix forming portion 24 reads the sets 
Q c and Q F stored in the element storage portion 23, 
discretization of the sets takes place using their elements 
such that a finite element method is applied to a Maxwell 
equation, and a matrix and column vector in simultaneous 
linear equations are formed (S12). Here, the Maxwell equation 
as an object is represented by the following expressions: 

Expressions 1 and 2 



Where E: electric field, H: magnetic field, e: dielectric 



VxH = (a + jo)i)E 
V x E = - jcojjH 



■••(1) 
-(2) 
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constant, |A: permeability, a: conductivity, co: angular 

frequency, V: differential operator (d/dx, d/dy, d/dz), X: 

outer product, j: imaginary number 

From these expressions, the following differential 

equation is obtained: 

Expression 3 
Vx pVx u-(Q 2 qu = 0 --(3) 

u represents an electric field, magnetic field, etc., and 
p and q represent characteristic values of dielectric constant 
and permeability. For example, when Expression (3) is 
obtained by substituting Expression (2) into Expression (1) 
and eliminating H, u represents an electric field, p represent 
the reciprocal of a permeability, and q represents a 
dielectric constant . 

When a finite element method is applied to an element V k 
formed in the former process described above, the following 
expression is obtained: 

Expression 4 
I (V x N ki )p(Vxu)dV-co 2 N ki qudV = {N k , xpWxu)-ndS 

-(4) 

Here, dV k represents the surface of V k and n represents an 
outward normal unit vector. N ki represents an interpolation 
function (also called a shape function, insertion function, or 
basis function) . When tetrahedral elements and linear 
equations are used, the interpolation function is given by the 
following expression : 
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Expression 5 



N kl = 

K 1 


1 1 | V Z J 




Nk7 = 


1 Z [ V J I 






|/ 3 |(AVi 2 - 


-AVI,) 


N k4 = 


I '4 |(W 




N k5 = 


|/ 5 |feVL 4 


-4.VI2) 


N k6 = 




-4VI3) 



Here, Li to L4 represent a volume coordinate system and 
I li I to I lei represent the length of a side Ij.. Moreover, 
regarding the detail of the volume coordinate system, 
Reference Document 2 (T. Itoh, G. Pelosi and P. P. Silvester, 
"Finite Element Software for Microwave Engineering," John 
Wiley & Sons, pp. 101-125, 1996) may be referred to. 

Fig. 6 shows components of an electromagnetic field u that 
are in contact with the sides of a tetrahedron. The 
components of the electromagnetic field u are arranged so as 
to be disposed along the six sides, and are represented by u x 
to u 6 . The electromagnetic field at point x inside the element 
V k is represented by the following expression using Expression 
(6) and the six components of an electromagnetic field: 

Expression 6 

«(*) = £ -(6) 

When Expression (6) is substituted into Expression (4) and 
applied to all of the elements of Q F , the following 
simultaneous linear equations are obtained: 

Expression 7 
A F u f = b F -(7) 

Here, u F represents vectors which are electromagnetic 
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fields disposed at the sides, and the dimensions are equal to 
the number of the sides as a general rule. The components of 
a matrix A F are formed by the following expression: 
Expression 8 

IIS ilCVxNJpCVxNJdV-rflNjNtdV] -(8) 

*=1 M j=l k k 

Furthermore, a vector b F is formed by the following 
expression . 

Expression 9 

tit LWvxpVxNJ-ndS ••■(» 

*=1 r=l j=) k 

Here, N F is the number of elements in Q F . 

The vector/matrix forming portion 24 also forms a matrix 
A c corresponding to Q c via the same procedure and preserves A F , 
b F , and A c in the vector /matrix storage portion 26. When u F is 
determined by calculating simultaneous linear equations in 
Expression (7), the electromagnetic field distribution can be 
calculated using Expression (6) to produce the solutions of 
Maxwell equations in Expressions (1) and (2). 

Next, the prolongation matrix forming portion 25 reads 
elements £2 F and Q c formed in the former process from the 
element storage portion 23 and, based on the geometric 
information, the prolongation matrix forming portion 25 forms 
a prolongation matrix P which relates the electromagnetic 
field vector u F in Q F to the electromagnetic field vector u c in 
Q. c (S13). When the prolongation matrix is used, as shown in 
the following expressions, one vector can be used to determine 
the other vector: 
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Expressions 1 0 and 1 1 
u F =Pu c -(10) 
U C =P t u F -(ID 

Here, t means a transposed matrix. 

Fig. 7 shows the relationship between an electromagnetic 
field vector u F in Q F and an electromagnetic field vector u c in 
Q c . The j component Uj C of the column vector u c is disposed at 
the side of an element of f2 c and the electromagnetic field u(x) 
at the position x in Q F can be represented using an 
interpolation function N c in Q c : 

Expression 1 2 

j 

The i component Ui F of u F is an average value of an 



electromagnetic field at the side li where the i component Ui F 
is disposed. That is, the i component ui F of u F is represented 
by the following expression: 
Expression 13 

«r-7nH*K* - (13) 

Here, ti represents a tangential unit vector of the side 
li. When Expression (12) is substituted for Expression (13) 
and compared with Expression (10), the ij component Pij of a 
prolongation matrix P is given by the following expression. 

Expression 14 

p ^\t\^ ixHdl "" (14) 

The prolongation matrix forming portion 25 forms a 
prolongation matrix using Expression (14) and preserves the 
prolongation matrix in the prolongation matrix storage portion 
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Next, the approximate value calculation portion 28 sets an 



electromagnetic field vector v using simultaneous linear 
equations according to a multi-grid method (S14). The initial 
value v F may be zero. Moreover, the details of the calculation 
of the simultaneous linear equations according to the multi- 
grid method may be referred to in Reference Document 3 (W. 
Briggs, V. Henson and S. McCormick, "A Multigrid Tutorial, " 
SIAM) . 

Next, the approximate value calculation portion 28 removes 
high-frequency components of error contained in v F by applying 
a stationary iteration method of simultaneous linear equations, 
such as the Jacob method, the Gauss-Seidel method or the SOR 
method, to Expression (7) (S15). Moreover, the details of the 
iteration method can be referred to in Reference Document 4 
(Satomi Hasegawa, Hidehiko Hasegawa, Kiyoji Fujino, "Iteration 
Method Templates", Asakura Shoten, 1996). 

Next, the approximate value calculation portion 28 maps a 
residual r F in Q F to a residual r c in Q c by using the following 
expressions in order to remove low-frequency components (S16). 

Expressions 15 and 16 



Next, the approximate value calculation portion 28 
determines solutions of simultaneous linear equations in the 
following expression (17) and calculates a correction vector e c 



approximate initial value v F to u F in order to calculate the 



r F = b F - 

r c =P'r F 




•••(15) 
•••(16) 



(S17) : 
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Expression 17 




■•■(17) 



In Expression (17), although a non-stationary iteration 
method, such as a direct method according to the Gauss 
elimination method, a CG method or a GMRES method is 
preferably used, since Expression (17) is a small matrix as 
compared to Expression (7) , the calculation can be performed 
in a short time. Moreover, the details of the iteration 
method can be referred to in the above-mentioned document 4. 

Next, the approximate solution correction portion 29 
improves the accuracy of the approximate solution by returning 
the correction vector of Q c to the correction vector of Q F by 
using the following Expressions (18) and (19) at (S18): 

Expression 18 and 19 



The approximate solution correction portion 29 makes high- 
frequency components attenuated by applying smoothing shown in 
step 15 a couple of times to the approximate solution (S19). 

Although the above processing in step S1 5 to S19 is the 
basic algorithm of a multi-grid method, since there are cases 
where the convergence is relatively slow, the convergence is 
accelerated by using another iteration method for simultaneous 
linear equations together as required (S20) . Specifically, 
the above-described processing is used as a pretreatment for a 
non-stationary iteration method, such as a CG method, a GMRES 
method or a GCR method, or the processing is accelerated by a 
residual cutting method, etc. Moreover, the details of the 
residual cutting method can be referred to in Reference 



e F = Pe c 
v F =v F +e F 



•••(18) 
•••(19) 
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Document 5 (Kikuchi et al . , JSME International Journal, 62-604, 
Series B p. 4076, December 1996). 

Next, the approximate solution correction portion 29 
evaluates the residual of Expression (15) in step (S21 ) . When 
the residual is large (S21 is no) , since the electromagnetic 
field vector does not converge, the process returns to step 
S1 5 and the following processing is repeated. Furthermore, 
when the residual has a sufficiently small value (S21 is yes), 
it is determined that the electromagnetic field vector 
converges, and the approximate solution v F is stored in the 
electromagnetic field vector storage portion 30 (S22) and the 
processing is finished. 

Fig. 8 shows one example of a rectangular waveguide as an 
analysis object of an electromagnetic field analyzer according 
to a preferred embodiment of the present invention. The 
dimensions of this rectangular waveguide are as shown in Fig. 
8 and the inside is filled with air. The electromagnetic 
field analysis is performed using tetrahedral primary elements. 

In Fig. 9, the calculation time to the number of elements 
when the rectangular waveguide shown in Fig. 8 is driven in 
TE1 0 mode of 2.45 GHz is compared between the case where a 
related direct method is used and the case where a method 
according to the present invention is used. As understood 
from Fig. 9, when the number of elements increases, the effect 
of the case where an electromagnetic field analyzer according 
to the present preferred embodiment is used is more pronounced 
and, when the number of elements is 206,266, the speed of the 
processing is about 17 times greater than that of the related 
direct method. 
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Fig. 10 shows one example of a patch antenna as an 
analysis object of the electromagnetic field analyzer 
according to a preferred embodiment of the present invention. 
The dimensions of the patch antenna are shown in Fig. 10, and 
the patch antenna is disposed on a dielectric substrate having 
a thickness of about 0.794 mm and a dielectric constant of 
about 2.2. The electromagnetic field analysis is performed by 
using primary tetrahedral elements. 

In Fig. 11, the calculation time to the number of elements 
when the patch antenna shown in Fig. 10 is driven at 7.0 GHz 
is compared between the case where a related direct method is 
used and the case where a method according to a preferred 
embodiment of the present invention is used. As is understood 
from Fig. 11, when the number of elements increases, the 
effect of the case where an electromagnetic field analyzer of 
the present preferred embodiment is used is pronounced and, 
when the number of elements is 88,445, the speed of the 
processing is about 1 5 times than that of the related direct 
method. 

As described above, according to the electromagnetic field 
analyzer of the present preferred embodiment, since a 
prolongation matrix P defining the electromagnetic field 
vector u F in Q. F related to the electromagnetic field vector u c 
in Q c is formed and an electromagnetic field vector is 
calculated by using this prolongation matrix, it is possible 
to perform an electromagnetic field analysis by using a multi- 
grid method using a non-nested mesh. 

Furthermore, since it is possible to perform an 
electromagnetic field analysis by using a multi-grid method 
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using a non-nested mesh, the time required for analysis is 
greatly reduced as compared to an electromagnetic field 
analysis using a direct method. 

While the present invention has been described with respect 
to preferred embodiments, it will be apparent to those skilled 
in the art that the disclosed invention may be modified in 
numerous ways and may assume many embodiments other than those 
specifically set out and described above. Accordingly, it is 
intended by the appended claims to cover all modifications of 
the present invention that fall within the true spirit and 
scope of the invention. 
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